Dimensional reduction in superconducting arrays and frustrated magnets 
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Some frustrated magnets and superconducting arrays possess unusual symmetries that cause the 
free energy or other physics of a D- dimensional quantum or classical problem to be that of a 
different problem in a reduced dimension d < D. Examples in two spatial dimensions include the 
square-lattice p + ip superconducting array, the Heisenberg antiferromagnet on the checkerboard 
lattice (studied by a combination of 1/S expansion and numerical transfer matrix), and the ring- 
exchange superconducting array. Physical consequences are discussed both for "weak" dimensional 
reduction, which appears only in the ground state degeneracy, and "strong" dimensional reduction, 
which applies throughout the phase diagram. The "strong" dimensional reduction cases have the 
full lattice symmetry and do not decouple into independent chains, but their phase diagrams, self- 
dualities, and correlation functions indicate a reduced effective dimensionality. We find a general 
phase diagram for quantum dimensional reduction models in two quantum dimensions with TV-fold 
anisotropy, and obtain the Kosterlitz-Thouless-like phase transition as a deconfinenient of dipoles 
of 3D solitons. 

PACS numbers: 74.50+r 74.72-h 



I. INTRODUCTION 

The subject of this paper is a set of problems in con- 
densed matter physics where unusual symmetries gen- 
erate three- or four-body interactions without two-body 
interactions, even though the underlying Coulomb inter- 
action is a two-body interaction. The examples that we 
study are drawn from models of frustrated Heisenberg 
magnets and complex superconducting arrays. Our fo- 
cus here is on a dramatic consequence of the four-body 
interaction that appears naturally in these systems: the 
apparent dimensionality of the system is different from 
its effective dimensionality, either for the system's phase 
transition or its ground state degeneracy, or both. Most 
possible experiments on the systems we discuss should 
measure correlations that are sensitive to the effective 
dimensionality. 

Some examples of this phenomenon arc relatively sim- 
ple and well understood. The 1/5 expansion approach 
to the 2D checkerboard Heisenberg antiferromagnet^ la- 
bels the infinite degenerate classical (5 — > oo) ground 
states by local Ising variables. Including 1/5 corrections, 
which come from the zero-point energy of spin- wave fluc- 
tuations, partially splits the extensive 5 — > oo degen- 
eracy and generates a four-spin energy function on the 
Ising variables. As reviewed in Section II, this 1/5 en- 
ergy function supports a ID degeneracy of ground states, 
even though the model is two-dimensional and symmetric 
x <-> y. However, as we confirm numerically, this effective 
model has a finite-temperature phase transition (an ex- 
ample of "order by disorder"—) as usual in 2D. Only the 
ground state degeneracy suggests a reduced dimension- 
ality, and we label this "weak" dimensional reduction. 

Most of the discussion in this paper is toward under- 
standing more subtle examples, both quantum and clas- 
sical, and in arbitrary dimension, where dimensional re- 



duction persists for all couplings. An example appearing 
in a 2D quantum model for p + ip superconducting ar- 
rays shows an exact self-duality identical to that of the 
2D classical Ising model, plus other signatures of 2D clas- 
sical or ID quantum physics. This model is one of the 
few cases where it is possible to obtain exact information 
on a two-dimensional quantum phase transition. 

There are several physical problems (Section II) where 
the above physics may be tested experimentally. The 
simplest experimental signature in lattice models is cor- 
relations that are highly directional: in general the long- 
wavelength correlations near a quantum critical point are 
approximately isotropic, but that is not the case for the 
critical points we discuss. Although the correlations have 
the full symmetry of the lattice (there is no decoupling 
into "chains" ) , they do not have the continuous rotational 
symmetry of ordinary critical points. 

For clarity, we give a classical, fully solvable example 
of dimensional reduction (arising in a model of supercon- 
ducting arrays), so that the basic idea and physical sig- 
nificance will be clear. This model has no finite-coupling 
phase transition; most examples in later sections have 
unusual classical or quantum phase transitions. The sim- 
plest example of what we mean by "dimensional reduc- 
tion" is the 2D classical theory on the square lattice with 
an Ising spin S{ = ±1 at each site and plaquette interac- 
tions: 
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where the sum is over elementary squares and the four 
spins are the corners of square □. In Cartesian coor- 
dinates 
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FIG. 1: A typical ground state configuration of Ising variables 
(open or filled circles) for plaquette interactions, as appear in 
a p ± ip superconducting array. Note that fixing one row 
and one column determines the rest of the lattice via the 
constraint of 0, 2, 4 up-spins per square. 



Section II reviews how this model gives the effective de- 
scription of p ± ip superconducting arrays^. 

This model has exactly the same free energy per 
site as the ID nearest-neighbor Ising model —j3E — 
KJ2i s i s i+i- To see this, note that just as in the ID 
Ising model all the bonds bi = s^s^+i can be specified 
independently (up to two spins at the boundary), in the 
2D plaquette model Q all face variables 
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can be specified independently (up to a single horizontal 
line and single vertical line of spins). 

Given arbitrary values of the spin on one horizontal 
line and one vertical line, there is a unique ground state 
of the system compatible with these values (Fig. [I}. Note 
that QJ differs from conventional Z2 gauge theory in that 
the Ising variables are defined on sites rather than bonds. 
There is a one-dimensional gauge group 5 , and the ground 
states are gauge-equivalent to the uniform state. This 
ground-state degeneracy may seem topological because 
the 2D model has a ID entropy, which scales with the 
boundary of the model. However, the ground-state de- 
generacy is relatively insensitive to the topology of the 
2D manifold, unlike the topological degeneracies in e.g. 
fractional quantum Hall systems. 

We instead use the term "dimensional reduction". 
This term was first used to describe situations where 
the free energy of one problem in dimension d is ex- 
actly equivalent to the free energy of another problem 
in dimension d + 2, as in the equivalence between the 
hard-sphere gas in dimension d and branched polymers 
in dimension D = (d + 2)£i For the 2D classical model 
(|TJ, the free energy per site is the same as that of a ID 
model, and there is a similar relationship of correlations. 
In the more exotic cases discussed in section II, there are 
exact properties like duality relations and transition lo- 
cations that arc equivalent for one model in dimension d 
to another model in dimension d + n, although the free 
energies are demonstrably inequivalent. 

The ordering that exists at T = in this classical 
model can be viewed as "bond order" , since all horizon- 



FIG. 2: The cross-like 4-bond correlation function. The lead- 
ing nontrivial term is the product of all the plaquettes above 
except the center plaquette (marked). 



tal bonds in the same column are either ferromagnetic or 
antiferromagnetic, as for all vertical bonds in the same 
row. Alternately, since every closed loop on the square 
lattice has an even number of bonds from each column 
and row, the order can be defined in terms of the Wilson 
loop (the product of Sj variables around a closed loop). 
Dimensional reduction occurs not just for thermodynam- 
ics but also for correlation functions: two vertical bonds 
on the same row have correlation function 



(6162) = (tannic)!* 1 -* 2 ! 



(4) 



just as for the ID Ising model. Two vertical bonds on 
different rows have zero correlation by symmetry. 

However, the four-point correlations do not factorize 
as one might expect. The correlation function of two 
horizontal bonds on the same column and two vertical 
bonds on the same row is (Fig. |2J 
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(tanhiql^-^l+l^ 4 !- 2 
(tanh^)- 2 (6 1 6 2 )( P 3P 4 ). (5) 



This lack of factorization indicates that the model does 
not simply decouple into a set of independent lower- 
dimensional systems. Experiments are typically most 
sensitive to the two-point correlations, which show a 
strong anisotropy that will be the most dramatic exper- 
imental signature of dimensional reduction. Usually lat- 
tice systems that respect lattice symmetry are spatially 
isotropic sufficiently close to a critical point, but this does 
not hold for the examples we discuss. 

Our focus in this paper is on physical problems that 
naturally generate plaquette interactions on the square 
lattice without ordinary two-point interactions. How- 
ever, there are materials of considerable current inter- 
est, like the high-temperature parent material La2Cu04, 
that combine four-spin and two-spin interactions on the 
square lattice. Neutron scattering measurements on this 
material 8 indicate that the ring-exchange coupling Jq is 
about 40 percent of the two-spin coupling J. While the 
precise dimensional reduction we study is specific to the 
absence of two-spin coupling, some indications survive for 
finite but small two-spin coupling; work in this direction 
for a bosonic model U(l) case is frA 
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A helpful way to understand the 2D quantum models 
we study is through their classical 3D analgoues: these 
are anisotropic combinations of two-point interactions in 
the imaginary time direction, with gauge-like four-point 
interactions in the spacelikc planes. Dualities in some 
of these anisotropic 3D classical models were studied in 
a series of papers in the 1980s by Savit and collabora- 
torsi2iiiii£. The self-duality in these models can be un- 
derstood roughly as follows, for the Z 2 case: since the 
Ising model in 3D is dual to Z 2 gauge theory, it is plausi- 
ble that a hybrid of Ising and gauge interactions should 
be dual to itself. Dimensional reduction in these models 
is related to the existence of an infinite but not exten- 
sive number of conservation laws. The solvable Kitaev 
models is an example of a 2D quantum model with an 
extensive number of conservation laws. The existence of 
an ID-infinite number of conservation laws is essential 
to some models of 2D fractionalized phases; the strong 
dimensional reduction models are ideal examples of this 
class for which exact nonperturbative results can be de- 
rived. 

The outline of this paper is as follows. Section II intro- 
duces three physical examples of partial or complete di- 
mensional reduction, and presents numerical results that 
confirm previous 1/5 results^ for the Heisenberg model 
on the checkerboard lattice, an example of weak dimen- 
sional reduction. Section III discusses dualities and phase 
diagrams of 2D quantum or 3D classical models that show 
strong dimensional reduction. In section IV we discuss 
the correlation functions of these theories in more de- 
tail, and in particular the sense in which the correlations 
are characteristic of a reduced effective dimensionality. 
Section V summarizes our main conclusions and some 
experimentally detectable signatures of dimensional re- 
duction in different condensed matter systems. We also 
discuss the relationship between these dimensional reduc- 
tion cases and the fully solvable Kitaev model. 



II. PHYSICAL REALIZATIONS 

This section discusses three different physical problems 
that lead to a reduction of effective dimensionality, cither 
for the ground state dimensionality ( "weak" dimensional 
reduction) or for both the ground state degeneracy and 
the quantum disordering phase transition ("strong" di- 
mensional reduction). The three problems are: ordering 
in superconducting arrays of Josephson-coupled p ± ip 
superconductors; the Heisenberg antiferromagnet on the 
checkerboard lattice; and the bosonic ring-exchange liq- 
uid introduced by Paramekanti, Balents, and Fisher 9 . 
The fundamental significance of the first two cases is 
that these problems generate four-spin interactions on 
the square lattice with zero two-spin interactions. The 
absence of two-spin interactions is essentially a conse- 
quence of symmetry in the first two examples, rather than 
of any fine-tuning. 

First we discuss frustrating phases in two-dimensional 



superconducting arrays made from grains of p ± ip su- 
perconductors like S^RuC^: the angular dependence of 
tunneling into such a grain leads to an effective Ising 
Hamiltonian with multiple-spin interactions. Such angu- 
lar dependence of tunneling provides the most powerful 
confirmation that the cuprate superconductors have d or- 
dei-li^. Complex order parameters like p±ip and d±id 
give new physics because the phase measured in Joseph- 
son tunneling winds continuously as a function of angle, 
while real order parameters only generate relative phases 
of or 7r. This angular-dependent phase causes several 
unusual collective effects in periodic two-dimensional ar- 
rays. 

It was shown in4 that, even in the absence of exter- 
nal magnetic field, there is a effective flux associated 
with each elementary loop that sometimes frustrates the 
Josephson energy. This flux is determined by the Ising- 
like time-reversal order parameters in the grains forming 
that loop. The Ising order parameter specifies which of 
two degenerate states related by time reversal exists in a 
particular grain: for example, d + id (a = +1) or d — id 
(a = —1). The classical Hamiltonian of such a system 
takes the Higgs gauge form 

H = -Ej ^2 cos(0i - - Aij(a-i,a-j)). (6) 

<ij> 

In the above 4>i is the phase of the superconducting order 
parameter in the ith grain and Oi is the Ising-like time 
reversal order parameter in that grain. The gauge field 
Aij for p ± ip superconductors is given by 

Aij = OiQi - aj6j, (7) 

where 9i is the angle between the crystalline y-axis of 
the ith grain and the normal of the interface from i 
to j. Hence the Josephson effect couples the discrete 
(Ising) variables <Tj to the continuous superconducting 
phase variables fa. 

In this paper we will concentrate on this model square- 
lattice of p ± ip grains for two reasons. First, its 
two-dimensional quantum phase transition is similar to 
the one- dimensional quantum Ising transition or two- 
dimensional classical Ising transition, and shows the same 
Kramers- Wannier self-duality^. Second, it shows a type 
of anisotropic "bond order" that also occurs in certain 
frustrated magnets. First consider the ground states of 
this problem. The frustrating phases disappear on the 
square lattice if, for every unit square, the four spins 
Si = ±1 satisfy 

7T 

— (si + S2 + S3 + S4) = 2im, neZ (8) 
which is equivalent to 

SlS 2 S 3 S4 = 1. (9) 

Hence the ground states of the system are exactly those 
of the gauge-like theory discussed in the introduction: 

0E=-K^s^ (10) 
□ 
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This differs from ordinary Z2 gauge theory in that the 
spin variables live on sites, rather than bonds. One re- 
sult is that the gauge group is one-dimensional rather 
than two-dimensional: the total number of ground states 
is 2 L * +L y for a system of size L x x L y . A typical ground 
state configuration is depicted in Fig. ^ Note that fixing 
the spins on one horizontal and one vertical line deter- 
mines all the remaining spins. Since each face can be cho- 
sen independently to have YlieD s i = the free energy 
per face (or per site) is — [3~ 1 \ogZ = — /31og(2coshif), 
just as for the ID Ising model, where all the bonds . 

We will focus on adding quantum effects to the sim- 
plified model QJ, even though this model is only exactly 
equivalent to the superconducting array problem for the 
ground states (an extended discussion of the classical 
physics of such superconducting arrays is ini). The justi- 
fication for treating the simplified model is as follows: one 
can imagine integrating out the superconducting phases 
to obtain an effective action for the Ising spins. At zero 
temperature, where the system is in a ground state, this 
still gives a nontrivial set of states, essentially on sym- 
metry grounds. Our treatment in this paper will start 
from the simpler model ; while the true spin action is 
presumably more complicated and includes long-ranged 
interactions between frustrated plaquettesA, it will have 
the same gauge symmetries as QJ, so a reasonable expec- 
tation is that similar physics will apply even in the much 
more complicated array model. 

The model becomes quantum-mechanical in the pres- 
ence of a transverse magnetic field: 

H = -K"£^»t (11) 

□ i 

Here again the K interaction is around a plaquette but 
now the spin is a quantum spin-half and the a are Pauli 
matrices. This model has an infinite number of conserved 
quantities, as the total S x spin on each row and column 
is fixed modulo 2»* In the superconducting array real- 
ization, the magnetic field h corresponds to tunneling 
between the two order parameters p ± ip, which will in 
the limit of strong tunneling induce a single real order 
parameter p x . Applying of pressure in Sr2Ru04 is found 
experimentally to drive the system toward a real p state, 
but the explanation of this effect is not yet known. This 
quantum Hamiltonian Ijlll) was discussed previously in 5 . 
There is an exact self-duality between the large K/h and 
small K/h phases, and a critical point at if = ft. In this 
paper, we will derive some new properties of this model 
and discuss how its physics is related to that of other 
examples of dimensional reduction. 

Another problem that naturally generates multiple- 
spin interactions without two-spin interactions is the 
1/5 expansion approach for highly frustrated magnets. 
Such magnets can have a large degeneracy of classical 
ground states, i.e., at 5 — 00, when the local spin is a 
classical vector. This degeneracy is split by quantum- 
mechanical corrections at order 1/5: for many cases, 
the effective Hamiltonian to order 1/5 is known to have 



multiple-spin interactions between Ising spin variables 
that parametrize the degenerate manifold of classical 
ground states. 

Here we will study one example of this physics in or- 
der to understand "weak" dimensional reduction. The 
2D checkerboard antiferromagnet (Fig. [3J , introduced as 
a qualitative model for 3D pyrochlore antiferromagnets, 
has an infinite classical ground state degeneracy. Al- 
though the checkerboard lattice is not known to describe 
any real 2D system, it is intended to model the corner- 
sharing geometry of the 3D pyrochlore lattice that de- 
scribes many important materials. These ground states 
can be labeled by Ising variables Sj = ±1 on a square lat- 
tice, subject to a constraint. The plaquettes with diago- 
nal interactions are like tetrahedra in the pyrochlore lat- 
tice, and around any such tetrahedron, the corner spins 
satisfy 

si + s 2 + s 3 + s 4 = 0. (12) 

Adding quantum fluctuations at order 1/5 splits the 
classical degeneracy, essentially because different classi- 
cal ground states have different spin-wave frequencies^ (a 
recent review of this approach is£) . The effective Hamil- 
tonian for these constraints istiSiii 

H = -—— y%iS2S 3 S4- (13) 

^ □ 

This correctly gives the ground states of the checker- 
board lattice problem, but is only an approximation for 
the excited states, since gapless excitations from non- 
collinear spin configurations are ignored. We study the 
Ising spin problem here to understand the role of the con- 
straint (|12fl . as the only differences between this problem 
and the previous model Q is the constraint. However, 
the ground-state physics of this checkerboard model is 
similar: there is a one-dimensional rather than a two- 
dimensional degeneracy of ground states^. Hence there 
is an effective reduction of dimensionality as far as the 
ground states are concerned. 

We have confirmed numerically that the finite- 
temperature physics of the two models is very differ- 
ent: even though the ground-state ordering is similar, the 
checkerboard model does have a finite-temperature phase 
transition^, unlike the superconducting array model. 
This conclusion is obtained by large-scale transfer-matrix 
calculations, as unlike the previous case we do not have 
an analytic solution for finite temperature. The key to 
the trivial solution of the previous case was that all the 
plaquette variables were essentially independent, like the 
bond variables in the ID Ising chain; this seems no longer 
to be the case once the constraint (|12f) is added. The 
transfer-matrix results on strips of increasing size (Fig.[3J) 
show a peak in the specific heat at kT « 0.95 J, and this 
peak sharpens with increasing strip size (Table I). 

Finally, we briefly review the physical picture intro- 
duced in* that leads to a bosonic ring-exchange model. 
The likelihood of strong ring-exchange terms in certain 
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FIG. 3: Geometry used for the transfer matrix calculations 
of Table I. Inclined strips of odd width are terminated by 
"tetrahedra" at both sides. 



Width 


Cpcak 


Tpcak 


c"(T pcak ) 

c 


AT 
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0.183 


0.922 


2.62012 


0.618 


11 


0.175 


0.950 


2.78854 


0.599 


13 


0.170 


0.960 


3.05586 


0.572 


15 


0.166 


0.957 


3.32221 


0.549 


17 


0.163 


0.945 


3.48748 


0.535 



TABLE I: Results of numerical transfer-matrix calculations: 
estimates for location and width of peak in specific heat for 
strips of different widths. 

two-dimensional systems has been emphasized in that pa- 
per: both spin models and superconducting array models 
can have many-body interaction terms of the same order 
as the usual two-body terms. If the two-body terms can 
be tuned to zero, as discussed in£, then the many-body 
problem can have unexpected physics, including an alge- 
braically correlated phase in 2+1 dimensions. This is a 
physical realization of the 3D Kosterlitz-Thouless phase 
discussed mi£. 

The Hamiltonian for this model on the square lattice 
is related to that for ordinary s-wave Josephson junction 
arrays with competition between Josephson and charging 
energies. However, instead of a two-grain term Ej = 
—J cos(4>i — 4>j), the Josephson coupling is postulated to 
couple the four spins around a fundamental square. The 
resulting quantum Hamiltonian for commensurate filling 
(an integer number no of bosons per site) is 

H = -J^cos(0° - 0° + 0° - +C^2(m - n ) 2 

□ i 

(14) 

Here 4>i and are conjugate variables forphase and num- 
ber of the bosons, and refer for i = 1,2,3,4 to the 
upper left, upper right, lower right, and lower left grains 
around a fundamental square. 

One point that the reader should keep in mind is that 
the same "strong" dimensional reduction underlies both 
this bosonic U(l) model and the Z 2 model (fill) . There is 
a family of Zjy models that interpolates between the two 
limits, and shows the same two-stage symmetry breaking 
for large N as the 2D U(l) model with AT- fold anisotropy. 
Some properties of the phase transitions can be deter- 
mined using the exact duality discussed in the next sec- 
tion. We find that even though two-point correlations 



in these theories are highly directional (which allows the 
physics we find to be measured), these theories do not 
factor into decoupled ID quantum problems: rather di- 
mensional reduction is an effective long-wavelength prop- 
erty of some but not all correlation functions. 

III. DUALITIES AND DIPOLE 
DECONFINEMENT TRANSITIONS IN 2D 
QUANTUM LATTICE MODELS 

This section starts by explaining how the 2D quantum 
models with four-site interactions introduced in the pre- 
vious section show generalizations of standard 2D classi- 
cal dualities. These models have the full symmetry of the 
2D square lattice. We show the duality first, for the Z 2 
case, directly in the quantum notation, then introduce 
the associated 3D classical models which were previously 
shown to be self-dual. We then derive a picture for the 
phase transitions of these 3D classical models in terms 
of a deconfinement of dipoles of solitons, which are dual 
to the bond variables that order across the transition. 
We find that unusual symmetries in the problem lead to 
a fixed point of lower dimensionality; these symmetries 
are analogous to the "sliding symmetry" that exists in 
a hydrodynamic model of quantum Hall striped phases. 
A recent preprint^ on this model finds that the sliding 
symmetry leads to an effectively ID fixed point. In this 
section, we will find a similar conclusion in models that 
start with the full symmetry of the square lattice. 

The self-duality is simplest in the quantum model Qllj l . 
Just as the duality of the 2D Ising model can also be done 
in the ID Ising chain with transverse fieldi^, this quan- 
tum duality is the same as that of the classical anisotropic 
model given below. First, let us define a dual lattice: i' 
is on the face surrounded by i, i + x, i + y and i + x + y. 
Define new variables on the dual sites by 

£>i' — a i cr i+x a 'i+y a i+x+ij- (15) 

Also we define 

sf,=n^- as) 

r<i 

By r < i we mean that both the x and y coordinates of 
r are not greater than those of i. 

The new variables S x and S z have the same algebraic 
relations as a x and a z , since 

(S x ) 2 = {S z ) 2 = 1, 

S XgZ + S Z S X = Q ( 17 ) 

We can also check the Hamiltonian in (|11|) can be rewrit- 
ten in terms of these new variables as 

H = —KSf, — h Sfi Sfi + %Sft + yS Z , +x+yy (18) 
i> □ 

which has the same form as the original Hamiltonian, 
except exchanging h and K . This is one example of the 
quantum version of the dualities. 
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FIG. 4: The classical anisotropic 3D problem that describes 
the quantum critical point of the model 1111 , The two types 
of interactions (shaded bonds) are plaquette interactions in 
the planes normal to z, and bond interactions along z. 



This quantum version confirms the order that devel- 
ops at the quantum phase transition* in the Z2 case. For 
example, when K/h > 1, the system is in an ordered 
phase with "bond order"—, when h/K > 1, the system is 
in a disordered phase, but it is in the ordered phase in 
the dual theory. The order parameter in the dual theory 
is the dual bond S*,S*, + ~, or S?,S?, +y . Written in origi- 
nal variables, these order parameters can be represented 

aS Hr(y)<i(y),r(x)=i(x) a r > or Ilr(x)<i(x),r(y)=i(y) a r • r ( x ) 

and r(y) are the x and y components of coordinate r. 
Non- vanishing expectation value of this order parameter 
means long range correlation of <j x . 

We now introduce a set of classical models in order to 
find a field theory description and phase diagram. The 
relation between quantum and classical models for the 
above Z2 case is given in detail in*; here we will simply 
introduce the classical 3D Z^ Villain models that were 
previously shown to be self-dual*V A slightly different 
self-dual N = 2 model (Fig. was defined in* to make 
the notation consistent with the Kramers- Wannier dual- 
ity of the 2D fsing model: for Ising variables s» on the 
cubic lattice, 



W 



CXp(-?C C y ^ ^i^i^x^i^y^i-^x-^-y ^ ^i^i-\-z 

□ 



(19) 

the duality relation smh(2K xy ) sinh(2isr z ) = 1 is exactly 
the same, but a high-temperature expansion shows that 
the free energy is not. These classical models have an 
asymmetry between one direction t and the other two 
(x, y): define the Z^v Villain model for a cubic lattice 
(x,y,t) as 



Qi+x+y — Qi+x — Qi+y + Qi 



(21) 



and the integer j xy -i is defined on the plaquette sur- 
rounded by q i+ x+y, q%+x, qi+y, qi- 

The coupling in the dual theory is related to that in 
the original theory by 



dual P — 



N 2 



(22) 



For the Villain model defined above, the dual theory has 
exactly the same form as the original theory; demonstrat- 
ing this self-duality is a simple if lengthy calculation^*. 
If the theory has two critical points, they are connected 
by the relations 



PW2 = N 2 /(2ttY 



(23) 



which is used in the following several sections. If it has 
one critical point, it must lie at the self-dual coupling 
j3 c = N/ (2ir). For higher dimensions*, we can just replace 
A xy by A xyz , the alternating sum of the eight corners of 
a fundamental cube, and most of the following can be 
generalized to higher dimensions. 

The anisotropic Z^r Potts model, which is a more di- 
rect generalization of our Z2 model£, has similar duality 
relations. The Zat Potts model is 



JV-l 

E 

g=0 



27f 27T 

Z = > ' exp[^ 22 cos( — A t(?i ) + @22 cos( — A xy qi)}. 



(24) 

And following the formulation in2* of the self-duality of 
the 2D Potts model, one finds that for TV = 2,3,4, the 
models are self-dual; but for N > 5, this series of models 
are not exactly self-dual. 

Now we turn to models with U(l) symmetry. In the 
pure cosine form l|14|) . there is not an exact self-duality. 
The reason is that the self-duality of Zjv model is a du- 
ality between the vertex operator and the soliton opera- 
tor; the vertex operator is contained in the U(l) model 
by adding to it some additional terms, which will make 
the U(l) model reduce to a limit of Z^ ones. And in a 
pure cosine U(l) model, there is no vertex operator in 
the Lagrangian, so it lacks a natural field that is dual to 
the solitons. The partition function of the U(l) Villain 
model is 

/oc o 
DO Y, exp[-|^(A t ^-27r^ ;l ) 2 



{j M =-oo} 



-^(A^-2^^) 2 ], (25) 



00 +00 „ 



{g=-Oo} {i M =-Oo} 



This anisotropic 3D model has a Kosterlitz-Thouless-like 
transition*** and the duality relation 0ip2 = 1- The 
rest of this section discusses this phase transition as a 
confinement-deconfinement transition of dipoles of soli- 
tons and obtains a phase diagram for the Zjv models. 
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Solitons (a.k.a. monopoles) and other topological de- 
fects have behavior that is quite sensitive to dimensional- 
ity. Since the previous section showed that there is a deep 
connection between quantum ring-exchange models in 2D 
and classical 2D models, it is natural to ask about how 
solitons are different from ordinary 2D quantum models. 
In the 2D classical XY model, the soliton (vortex) can 
cause a phase transition, as can be understood from the 
famous argument of Kosterlitz and Thouless. A single 
vortex of integer winding number n has energy 



ft 



d 2 k — 

■n/L k 

7rn 2 /31n(L/a). 



(26) 



Here 1/a and 1/L are ultraviolet and infrared cut-offs. 
For simplicity consider the case with \n\ = 1. The vortex 
can take any place within the 2D plane, so it has roughly 
(L/a) 2 different positions. The entropy of a single vortex 
is thus about 21n(L/a). So the free energy goes like 



F = (n/3 -2)]n(L/a). 



(27) 



When j3 < ft c = 2/ir, the free energy will favor deconfmed 
vortices. 

This phase transition does not occur in isotropic 3D 
models. Because the energy of a Dirac monopole goes 
as fd 3 kl/k 2 , it is convergent in the infrared limit, 
which means the entropy will always dominate; there 
is no phase transition and the system is always in a 
monopole deconfmed phase. We have seen that certain 
3D anisotropic classical models have many properties 
that are analogous to a lower-dimensional model. Hence 
it is not too surprising that there can be a Kosterlitz- 
Thouless-like phase transition in the U(l) case: the 
existence of such a transition was first discussed for 
anisotropic 3D classical models Its appearance in 

2D quantum models and its physical consequences were 
studied irA Our goal in the following is to develop a 
physical interpretation of this transition similar to the 
vortex-unbinding picture of the 2D Kosterlitz-Thouless 
transition. To our knowledge, this discussion in terms of 
dipoles of solitons is new. This picture will then be used 
to obtain the phase diagram of this model with TV-fold 
anisotropy. 

The energy of a single soliton in the U(l) case, calcu- 
lated using the Gaussian part of the U(l) theory (|25|l . 
is 



2^ 



dk. 



1 



2tt k 2 k 2 + k 2 



4 Jir/L > 



/L 



7T ky 



/3(m-) 2 . 
a 



(28) 



The (2ir) 2 in front of this equation is inserted because 
that the charge of soliton is 2ir. The entropy is modified 
to S = 31n(L/a). Hence the free energy is 



F = /?(ln-^) 2 



31n(L/a) 



(29) 



In the limit of large L, the energy always dominates, and 
we can never get a deconfmed soliton. 

However, there is another way to drive a phase tran- 
sition. Instead of considering the single soliton, let us 
consider a dipole of solitons. By dipole we mean that mi 
takes values n and —n on two sites connected by a bond. 
For example, nij = 1, and mi +y = — 1. The energy of a 
fundamental dipole n = ±1 is 



(2 " } 2 J ^ J ^ 



dk 7 



k 2 



27T k 2 k 2 y + k 2 



tt/L 



27r) 2 ^lnL.(30) 



ir/L 



7r k x ky 



In the formula above, we have taken a = 1; because ev- 
ery derivative contains a in front of it, so if a =/= 1, we 
can redefine the lattice in every direction, x' — ax, and 
the theory should be the same with the case with a = 1. 
Now we can see that the energy of dipole is not infrared 
convergent, it also has a logarithmic divergence, of the 
same order of entropy. This indicates that a Kosterlitz- 
Thouless-like transition is possible, but there is an im- 
portant difference in that the energy calculation is con- 
vergent in the y direction, and hence cutoff-dependent. 

The entropy in this case is no longer 3 ln(L/a), because 
we saw that the energy is convergent in y direction, so 
the location of the dipole in y direction does not mat- 
ter. Actually, as the y-directed dipoles are uncorrelated 
if their y coordinates differ, we can say that their entropy 
comes only from their location in the x-z plane, so the 
entropy is of order 2 ln(L/a). Now we will find that there 
is a critical value of (3 beyond which the dipoles are con- 
fined to be a quadrupole; otherwise there is a deconfmed 
dipole phase. The procedure for obtaining this value is 
more subtle than in the 2D case. 

Dipoles with larger moment have higher energy, so they 
are more easily confined; this applies both to dipoles with 
higher charge and to dipoles which have bigger distance 
between their positive and negative charges. We can see 
that in the calculation above we have put in a cut-off of 
7r for k y . However, this only occurs because this model is 
so anisotropic that the y-directed dipole only has corre- 
lations in x and z directions: two y-directed dipoles can 
only be correlated when they are at the same x-z plane. 

With this in mind, we want to find a effective field 
theory for (|25|l that can describe this model near the 
new Kosterlitz-Thouless critical point. In the 2D case, 
in order to describe the vortices, we can replace the vor- 
tex coupling in the dual theory i2inn4> by cos(27T -m(f)) . 
When we expand the cosine term order by order, we 
get the vortex gas with charge 2um. And because the 
vortex with charge 2~k is the one with lowest scaling di- 
mension, we need only consider the \m\ = 1 case, be- 
cause near the critical point, vortices with higher charges 
are all irrelevant. In our 3D model, the single soli- 
ton is never relevant; we only want to find a term 
that can reproduce the dipole gas. The most relevant 
dipoles consist of two charges ±1 located on neighboring 
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sites. So we should replace the soliton-field coupling by 
cos(2nd x <j)) + cos(27T(9 y 0). 

Hence we consider as a trial theory of this critical point 
the Lagrangian density for dipoles 

+a cos(2nd x (j)) + a cos(27r<9 a 0). (31) 

When the two cosine terms are relevant, the dipoles de- 
confine, just as the vortices deconfine in the 2D case: 
just as cos(27r0) generates vortices in the usual 2D case, 
cos(27r92;0) generates an x-directed dipole. Our approach 
will be to argue that because of a gauge symmetry in this 
Lagrangian, the effective theory is not 3D but 2D, and 
then calculate various properties including the phase di- 
agram from the 2D theory. 

Let us see when the cosine terms are relevant for 
dipoles in the x direction. We want to compute the corre- 
lation function (exp(2iird x (f>(ri)) exp(— 2iitd x (f)(r2))) ■ For 
simplicity, we assume that z% = z%. This is just a simple 
Gaussian integral, and easy to evaluate as long as we pay 
attention to the cancelling of the infrared cut-off. This 
requirement in the 2D case gives rise to the result that 
the total charge in any order in perturbation theory is 0, 
so the total charge is conserved to be 0. One finds that 
cut-off independence implies that the correlation can be 
non-zero only if x% — xi'- dipoles are uncorrelated along 
the direction of the dipole moment. 

This arises from a ID gauge symmetry in the field 
theory fJTJl. Without the cosine terms, the Lagrangian 
is invariant under transformations 0-^0 + fi(x) + 
fiiy), which means there is a particle number conser- 
vation along each row and column. If we identify <p 
as a phase, then exp(2?7T(/>) is just the particle creation 
operator. So any two-point correlation function like 
(cxp(2i7r0(ri)) exp(— 2iir<f>{r2))) is always zero because it 
breaks the conservation along some columns or rows; this 
is another way to see that the single soliton is always ir- 
relevant. The first nonzero correlation function is 

(exp(2i7T0(xi, yi)) exp(-2i7r0(a;2, yi)) 
x e-xp(2iir<j)(x2, yi)) exp(-2i7r^(xi, y 2 ))}- (32) 

This is a four-point correlation function and the conser- 
vation requires that the four points are four vertices of a 
rectangle. If we take xi = x± + a, we come back to the 
case of dipole correlation function, and the x coordinate 
of two dipoles must be the same. 

The field theory in (|31|l seems anisotropic, but when 
calculating the correlation function of two x-directed 
dipoles (which must have the same value of y), one gets 
isotropic behavior in the z and y directions. A naive cal- 
culation of the correlation function for any y\ — y2 and 
Z\ — Z2, as in the 2D case, gives 

~( 2 A 2 2 )^+-- (33) 
y + v z 



This seems to predict that the phase transition occurs 
when the scaling dimension of the dipole operator is 2, 
at 

Bl_ =2 ^^ = - (34) 

However, this calculation is incorrect: when we calcu- 
lated the scaling dimension of the dipoles, we were us- 
ing the Gaussian theory, which means we neglected the 
quadrupole fluctuations on microscopic scales. In the 2D 
XY model, one can neglect the dipole fluctuations be- 
cause the scaling dimension calculation is insensitive to 
microscopic details, but here the cutoff dependence of the 
final answer indicates that microscopic fluctuations will 
modify the critical coupling /3|. 

We have not been able to calculate all the effects of 
quadrupole fluctuation directly. Noticing that the corre- 
lation functions of dipoles are only nonvanishing in one 
row or column, depending on the direction of the dipole, 
and also that the correlation is isotropic in either the z-x 
directions or z-y directions, we guess that the field the- 
ory for the dipolar transition is actually the 2D classical 
U(l) theory. This conclusion would be precise if all cor- 
relations were confined to these planes, as a field theory is 
determined by the set of all its correlations; instead there 
are higher-order correlations that are not purely planar, 
as in the Introduction, but these seem to be irrelevant at 
the dipolar transition. 

Under this hypothesis, the leading correlations of the 
anisotropic 3D theory are described by a set of decoupled 
planar effective Lagrangians for dipoles: 

C = Yd 5>'^ 2 - Yp E( 9 ^) 2 + «cos(2^). (35) 

Here cos 2-Ktp generates y-directed dipoles that drive the 
transition when this term is relevant, and both sums are 
over one plane normal to y. There is a similar Lagrangian 
for x-directed dipoles, of course. The renormalized cou- 
pling B determines the planar scaling dimension of the 
dipole operator as in the 2D case: 

A 2 = 7T/3. (36) 

For > f3 c = 2/tt, there is a Gaussian phase with effec- 
tive Lagrangian (now for y-directed bonds) 

£ = §E(^) 2 -iE(^) 2 - ( 37 ) 

This hypothesis for the field theory immediately sug- 
gests a phase diagram for the self-dual Zjv models intro- 
duced above; since the effective field theory for dipoles is 
a 2D classical theory, the phase diagram is very similar 
to the classical 2D case. First recall the effect of TV-fold 
anisotropy on the 2D XY model. The partition function 
of the 2D Zjv model is 

Z= Y, e X p[/?5>os(— A„ ft )]. (38) 

{q=0} i 
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From now on, we assume that the lattice constant in all 
the directions is 1. Following custom, we replace the 
models with a simpler model from the same universality 
class, the Zjv Villain model: 



Z = 



N-i +00 „ 2tt 

E E cx p[-fE(^F 

{9=0} {i M =-°o} i 

(39) 

The Zjv model is self-dual, with the duality relation^ 



0102 



4^' 



(40) 



Renormalization group analysis predicts that, when 
N > 4, this model has two critical pointsSi, and the crit- 
ical (3 C of the two critical points have the duality relation 
0\0^ = N 2 /(Att 2 ). One of the two critical points, say 
corresponds to the change in relevance of the symmetry 
breaking term, and the other one (3% can be viewed as 
the Kosterlitz-Thouless phase transition, with 

01 = -- (41) 

At low temperature, vortices are confined, the chemical 
potential of vortices flows to under the renormaliza- 
tion group. But at high enough temperature, (3 < (3^, 
vortices become relevant. So the self-duality is also a du- 
ality between vertex and vortex operator. The (27r) 2 in 
the denominator of (|4*U|> is due to the fact that the charge 
of a vortex is 2tt in these units. 

The cases N < 4 have only one critical point, whose lo- 
cation is fixed by the duality relation, (3 C — N/(2tt). The 
reason is that when we calculated the scaling dimensions 
of cosine terms, we neglected the contribution of vortices; 
this is only consistent when the phase transition caused 
by the cosine term takes place in a vortex-confined phase, 
i.e. f3f > @2- This only happens when N > 4, when we 
have a Gaussian phase. In the case of N < 4, the temper- 
ature of the Kosterlitz-Thouless point is lower than the 
temperature that make the cosine term irrelevant, i.e. 
0i < 02i which means that there is no Gaussian phase. 

Now carrying over these results to our problem is very 
simple: instead of spins we have bonds, and instead of 
vortices we have dipoles. First, the phase transition due 
to dipoles take place when the scaling dimension of the 
dipole operator becomes 2, from the last formula in last 
section we know that this happens at the renormalized 
coupling 

C 2 = (42) 

We can calculate the scaling dimension of the bonds in 
terms of the Gaussian theory, and get the result 



Ax 



N 2 
Att(3" 



The phase transition happens when 

8tt' 



01 



(43) 



(44) 
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FIG. 5: The phase diagram of the 3D anisotropic classical 
Zjv models: when N < 4 there is one critical point, and when 
N > 4 there are two. 



If (3 1 — 02i we S e t the critical value of N, which gives 
us N c = 4, the same result as in the classical 2D case. 
This is consistent with the phase diagram obtained 2 ^ by 
numerical simulation. 

The bond-dipole duality is also applicable to another 
kind of model. There are numerical resultsii on 3D clas- 
sical plaquette Zjv models of the form 



H = -K[cos( 



2wA xy qi 27rA yz qi 



N 



-) + cos( 



+ cos( 



N 

2-K/\ xz qj 
N 



)}■ (45) 



The bond-dipole picture is also applicable here: this 
model also has a similar gauge symmetry and the or- 
dered phase has a similar "bond order"—, but now with 
bonds in three different directions in the classical model, 
rather than bond order in the x and y directions and or- 
dinary (spin) order in the z direction. We therefore guess 
that the effective Lagrangian of bonds is the same as be- 
fore but with also z-directed dipoles propagating in the 
x-y planes. From this Lagrangian and its dual form, the 
critical N c is still 4, consistent with numericsii. Hence 
the phase diagram JSJ is fairly universal for such models. 



IV. CORRELATION FUNCTIONS 

The experimentally accessible physics of the above 
quantum 2D models is contained in their correlation func- 
tions, which exhibit a strong d-wave-like anisotropy even 
near the critical point. The very simple model presented 
in the introduction is a good example of how the cor- 
relation functions behave in the disordered phases: the 
two-point correlations of bonds are localized along rows 
or columns because of the gauge-like symmetry, but the 
four-point correlations do not factorize trivially. In the 
ordered phases, there is spontaneous symmetry breaking 
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and hence the correlations need not be localized; there 
are also Gaussian phases intermediate between order and 
disorder. 

Sections I and III have shown that the conservation 
laws or symmetries require that the 2-point correlation 
function between two vertical bonds, say, should vanish 
when the 2 bonds are not along the same row, which 
was crucial in justifying the lower-dimensional effective 
theory. However, the 4-bond correlation functions are 
nontrivial even in the disordered phase, just as in Section 
I. Taking the 1 2 case as an example, the classical 3D 
partition function can be written as 5 

Z = 2J exp(K 2J SiSi+zSi+fiSi+s+Q + H 2J SiS i+i ). 

{s} □ b 

(46) 

First let us calculate the equal-time 2-bond corre- 
lation functions (s Xl , y s Xlt y + iS X2ty s X2t y + i) in a high- 
temperature expansion. The first term is 

G (x u x 2 ) = ifl* 1 -^! = e-l^-^l/S. (47) 

Here £ = — 1/ (In K) . This leading term comes from ex- 
panding the Hamiltonian along the shortest path con- 
necting the two bonds. The next leading order comes 
from adding one bend to this path, and gives 

G 1 =2K 4 \x 1 -x 2 \G (x 1 -x 2 ). (48) 

Now let us consider one type of four-bond equal-time 
correlation: 

G 4 (x 1 ,x 2 ,y,x,y 1 ,y 2 ) = {b xuy b X2 ). (49) 

Here p denotes a horizontal bond, and b a vertical bond. 
Provided that the shortest paths intersect, as in Fig. [21 
the leading term in the high-temperature expansion is 

G 4 (x 1 ,x 2 ,y,x,y 1 ,y 2 ) = —^G (xi, x 2 )G (yi, y 2 ) + ■■■ 

So the non-factorization is already present in the leading 
term. We also checked the next term: 

G i (x 1 ,x 2 ,y,x,y 1 ,y 2 ) = ■^G (x 1 ,x 2 )G (y 1 ,y 2 ) 

+10K 4 G (x 1 ,x 2 )G (yi,y 2 ). (51) 

Hence the 4-bond correlation function is not simply 
l/K 2 xG(xi,x 2 )G(yi,y 2 ). The physical picture in the Z 2 
case^ is that the bond in the disordered phase is the low- 
est mobile excitation, but normally the bond can only 
move transverse to its direction (i.e., horizontal bonds 
hop vertically). However, when two vertical bonds mov- 
ing in the same row has opposite momentum, they can 
interact and scatter into two horizontal bonds; the hori- 
zontal bonds can scatter with each other and become two 
vertical bonds. Then in the Hamiltonian there is an ef- 
fective interaction between this two kinds of excitations, 
and this kind of excitation make the 4-bond correlation 



function nontrivial. Calculating the equal-time correla- 
tion of four horizontal bonds, say, gives the product of 
two-bond correlations to leading order, but the next term 
does not factorize. 

The correlations in the Gaussian phase and their con- 
nection to conductivity were discussed irA The correla- 
tions in the long-range ordered phase are quite different: 
in this phase, the expectation value of the bond operator 
is nonzero, which means the "gauge" symmetry is broken. 
There is much less restriction on the physical correlation 
functions. For example, {b Xl yi b Xa ya ) can be nonzero for 
different y\ and y 2 . The physical reason for this is that in 
the ordered state is a condensate of the bond variables, 
and can be regarded as a reservoir of bonds. So the two 
vertical bonds in different rows can interact with each 
other by exchanging bonds with the condensate, and the 
two-bond correlation function is highly sensitive to the 
transition between disordered and ordered phases. 

V. CONCLUSIONS AND OTHER 2D 
QUANTUM MODELS 

We have shown that frustrated magnets and super- 
conducting arrays naturally generate plaquette interac- 
tions on the 2D square lattice, and that such interac- 
tions effectively reduce the system's dimensionality, even 
though the models have the full symmetry of 2D square 
lattice. This dimensional reduction can be an exact 
equality of free energies or an exact equivalence of du- 
ality relations, or ("weak" dimensional reduction) just 
a ground-state/infinite-coupling property, depending on 
microscopic details. The correlations in these theories 
exhibit the reduced dimensionality, and these theories 
provide a way to promote phenomena from one quantum 
dimension, like exact dualities and fractionalized phases, 
to higher dimensions. 

There is a somewhat similar model in two quantum 
dimensions that is exactly solvable^: for spin-half vari- 
ables a now defined on bonds of the square lattice, 

H = -K 1 J2 o-\e z 2 olo\ -K 2 J2 °l°t°l°t- (52) 
□ + 

The first term contains the four bonds around a plaque- 
tte, while the second contains the four bonds meeting 
at a vertex. These two terms commute with each other, 
unlike the two terms in our Z2 model (|11|) . and as a 
consequence the model ()">2t has a 2D infinite set of local 
conserved quantities and is exactly solvable for all K\ and 
K 2 . The models studied in this paper can be thought oS 
as intermediate between the generic 2D quantum model 
with no local conservation laws, and the Kitaev model 
with an extensive set of local conservation laws: "strong" 
dimensional reduction seems to depend on the existence 
of enough conservation laws that the model's effective 
dimensionality is reduced, but not to zero. It is quite 
remarkable that this exotic behavior occurs in models 
motivated by experiment. 
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